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Abstract 

Q H | We study the incorporation of QCD effects in the basic electroweak corrections 

Af , Ar w , and Ar. They include perturbative 0(aa s ) contributions and ti threshold 
effects. The latter are studied in the resonance and Green-function approaches, in 
the framework of dispersion relations that automatically satisfy relevant Ward iden- 
tities. Refinements in the treatment of the electroweak corrections, in both the MS 
and the on-shell schemes of renormalization, are introduced, including the decoupling 
of the top quark in certain amplitudes, its effect on e 2 (m z ) and sm 2 9 w (m z ), the in- 
corporation of recent results on the leading irreducible O(o 2 ) corrections, and simple 
expressions for the residual, i.e. "non-electromagnetic", parts of Af, Ar w , and Ar. 
The results are used to obtain accurate values for m w and sin 2 '6 w (m z ), as functions 
of mt and m H . The higher-order effects induce shifts in these parameters comparable 
to the expected experimental accuracy, and they increase the prediction for mt de- 
rived from current measurements. The MS and the on-shell calculations of Ar, in a 
recently proposed formulation, are compared and found to be in excellent agreement 
over the wide ranges 60 GeV < m H < 1 TeV, m z < mt < 250 GeV. 
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1 Introduction 



It has been recently shown that the Ss method of renormalization provides a very con- 
venient framework to incorporate higher-order corrections to Ar [fj] arising from resumma- 
tion of one-loop effects. These include not only leading logarithms of 0((a ln(m z /mj)) n ), 
where m/ is a generic fermion mass, and subleading logarithms of 0(a 2 ln(m z /mf)) [[J, 
but also terms of 0((am 2 /m 2 v ) n ) |], f|, |J . The reason can be traced to the fact that in 
this method one essentially subtracts the divergent parts of the amplitudes. In contrast 
with other approaches, this procedure circumvents the introduction of mass singularities 
and 0(am 2 /m 2 v ) terms via the finite parts of the counterterms. As a consequence, the 
renormalized perturbative expansion has a structure very similar to that of the bare theory, 
where resummation of one-loop effects is easy to implement [I], |J]. There are, of course, 
irreducible two-loop contributions of 0((am;/m^) 2 ). As discussed in Refs. these 
can be gleaned from Ref. || and the more recent work of Ref. on the 0(a 2 ) corrections 
to the p parameter. 

It has also been recently shown || that it is possible to derive a simple and accurate ex- 
pression for Ar, within the on-shell method of renormalization ||, which contains the same 
leading and subleading contributions described above. On the other hand, the irreducible 
two-loop contributions of 0(a 2 m 2 /m^) have not been computed, so that both the Ws and 
on-shell calculations of Ar become uncertain at this level of accuracy. In fact, one can see 
that the difference between the two calculations and their inherent theoretical uncertainty 
due to the neglect of higher-order electroweak corrections start with subleading terms of 
0((a/TTS 2 )(c 2 1 s 2 )x t ) , where x t is the leading correction to the p parameter (cf. Eq. (|17b| )). 
or, equivalently, of C((c 2 /s 2 )(a 2 /27r)(fi: t /27r)), where a 2 = g 2 /^ and K t = (G t ) 2 /47r are 
the SU(2) and the Yukawa couplings of the top quark, respectively. These are very small 
for low m t and are expected to be ~ 8 x 10 -4 for m t = 250 GeV. (As we will see, over the 
large ranges m z < m t < 250 GeV, 60 GeV < m H < 1 TeV, the actual numerical evaluation 
of the on-shell expression of Ref. || and the ms calculations show a very small difference, 
reaching a maximum value of 2.5 x 10~ 4 at m t = 250 GeV and m H = 1 TeV, a very precise 
agreement which may be somewhat fortuitous.) This uncertainty is to be compared with 
an estimated error of ±9 x 10 -4 originating in the 0(a) contributions of the first five quark 
flavors |, [10 . 



In order to set the stage for our discussion, it is convenient to recall at this point some 
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of the basic relations of the on-shell and ms frameworks [|I|^,|TT|,|T2"||: 

A 2 

m 2 w (l - Ar)' 
A 2 

m 2 w {l-Af w y 
A 2 

m|(l — Ar) ' 
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where m w and m z are the physical masses, A = (7ra/(v / 2G At )j"" = (37.2802±0.0003) GeV, 
s 2 , s 2 , and c 2 are abbreviations for sin 2 ^ = 1 — m\ f jm\ and the ms parameters sin 2 9 w (m z ) 
and cos 2 # w (m z ), respectively, and Ar, Af w , and Ar are radiative corrections. The 't Hooft 
mass, /i, has been set equal to m z in Eqs. (|]) and (||). 

It follows from the analysis of Ref. |TB| that m z and m w in Eqs. ([!]) and (|2|) can 
be identified, phenomenologically, with the masses measured in current experiments and, 
theoretically, with the definition m\ = m\ + T 2 , where s = m\ — im 2 T 2 is the relevant 
complex pole position ||14)| . The latter is given by s = itlq + A(s), where mo is the bare 
mass and A(s) the conventional self energy, which includes tadpoles, tadpole counterterms, 
and, in the Z° case, 7Z mixing effects that start in 0(a 2 ). On general grounds, it is 



expected that m 2 and T 2 , and therefore mi, are gauge-invariant to all orders [|T^,|T5[. Over 
a large class of gauges, including those in which explicit calculations have been carried 
out, mi differs from the "field-theoretic" or "on-shell" definition, m 2 = m 2 , + 7teA(m 2 ), by 
gauge-dependent terms of 0(a 3 ) |T3|| . Because contributions of this order are well beyond 
the accuracy that may be achieved in the foreseeable future, the replacement of m by the 
more rigorous definition mi does not require a modification of the radiative corrections Ar, 
Af w , and Ar. Using the expression 
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Ar) 
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equivalent to Eq. (|IJ), the analogous one with s 2 — > s 2 and Ar — > Ar, equivalent to Eq. (|3|), 
and the accurately known value m z = (91.187 ± 0.007) GeV [[U|, the corrections Ar and 



Af lead to precise evaluations of s 2 and s 2 , as functions of m t and m H . These, in turn, can 
be compared with other determinations of s 2 and s 2 to constrain the value of m< and, in 
the future, that of m H . They are also important input parameters in the prediction of the 
Z° partial widths and on-resonance asymmetries, as some of these observables depend very 
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sensitively on the weak mixing angle (see, for example, Ref. [|17| ). More generally, the basic 
corrections Ar, Af w , and Af are frequently employed to verify the Standard Model (SM) 
at the level of its quantum corrections and in searches for signals of new physics |T8|-pT| . It 
was also explained in Ref. [|l|] how Af and Af w , relevant corrections in the Is framework, 
can be employed to evaluate the on-shell quantity Ar. 

The aim of the present paper is to incorporate the leading QCD corrections in the calcu- 
lations of the radiative corrections Ar, Af, and Af w . We also introduce some refinements 
in our previous treatment of the electroweak corrections [IJ. 

The relevant QCD contributions occur in the vacuum-polarization functions associated 
with the and Z° bosons and have been extensively discussed in the literature p2|-|3"T 



In particular, the QCD corrections involving the (t, b) isodoublet are known to be significant 
for large rrif. There are actually two types of effects that may be classified as perturbative 
0(aa s ) and threshold contributions. In the literature, the latter are frequently referred to 
as "non-perturbative." The perturbative 0(aa s ) components have been studied with two 
different methods: i) direct evaluation of the two-loop Feynman diagrams in dimensional 
regularization, an approach that goes back to the pioneering work of Djouadi and Verzeg- 
nassi 152]; ii) calculation of the imaginary parts and computation of the full amplitude by 



means of suitably defined dispersion relations (DRs) p3|-|26||. It has been shown pO|,pT| 
that the two approaches are equivalent in the evaluation of the perturbative contributions 
to Ar and Ap, a welcome fact. On the other hand, the Is scheme is implemented in the 
framework of dimensional regularization and, for a full determination of the subtraction 
constants, one must appeal to method (i). 

Threshold effects on the absorptive parts of the self energies have, in turn, been discussed 
in two different approaches: a) in Ref. |24| the contributions of densely spaced, narrow 
quarkonium resonances were taken into account on the basis of a specific quark- ant iquark 



potential; b) in Refs. [p7, 28[| one considers the imaginary part of the Green function for 
the non-relativistic Schrodinger equation that characterizes the tt system near threshold. 
The latter formulation effectively resums the contributions of soft-gluon exchanges in the 



ladder approximation (see also Ref. p9 



For sufficiently low mt there should be, near threshold, a rich spectrum of distinct non- 
relativistic states bound by strong long-range forces and the approach (a) is very natural. 
For increasing m t , however, the weak decay of a single top quark inside the bound states 
becomes important and, for m t ^ 130 GeV, the partial width of t — > W + b is so large that 
the revolution period of a tt bound state would exceed its lifetime. As a consequence, the 
individual resonances lose their distinctiveness and are smeared out to a coherent struc- 
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ture P7| ~ P5| . In that regime, the Green-function method is more appropriate. In summary, 
one expects approaches (a) and (b) to be preferable for lower (^130 GeV) and higher 
(>130 GeV) values of mt. Both formulations deal directly with the absorptive parts of 
the amplitudes. To obtain the real parts is then necessary to employ DRs. This was done 
in detail in Ref. |H] using DRs for the vacuum-polarization functions directly constructed 
from relevant Ward identities |[30|| . In conjunction with very plausible assumptions con- 
cerning the asymptotic behavior of the threshold effects for large q 2 , this procedure leads 
to specific results for the real parts. 

In Section 7, we compare four different calculations of m w and sin 2 9 w (m z ), as functions 
of m t and m H : i) only electroweak corrections; ii) electroweak plus perturbative 0(aa s ) 
corrections; Hi) the above, plus threshold effects in the resonance approach; iv) same as 
(Hi) with threshold effects in the Green-function framework. This allows us to demonstrate 
the magnitude of the QCD corrections and, at the same time, to separate the threshold 
effects from the more established perturbative 0(aa s ) contributions. Although the two 
approaches to treat the threshold effects are certainly not identical, we find the encouraging 
and fortunate result that their overall effects on m w and sin 2 9 w (m z ) are quite close over 
the entire range m z < m t < 250 GeV. 

As mentioned before, aside from incorporating the QCD effects in the relevant self ener- 
gies, we introduce some refinements in our treatment of the 0(a) electroweak corrections. 
Conceptually, the most interesting modification is a slight change in the definition of the 
fundamental couplings e 2 (m z ) and sin 2 9 w (m z ) of the Us scheme, which is introduced in 
order to make them essentially independent of heavy particles such as the top quark or 
unknown massive excitations. In the case of sm 2 9 w (m z ), we follow a convention recently 



proposed at the one- loop level by Marciano and Rosner |32|,[33l], and explain how to extend 
it when 0(aa s ) corrections are included. We emphasize that these modifications in the 
definitions of the fundamental Ws parameters do not affect, to the order of the calculations, 
the relations between physical observables because they are compensated by corresponding 
changes in the appropriate radiative corrections. A second change is that we use an up- 



dated calculation by Jegerlehner |T(J for the contribution of the first five quark flavors to the 
photon vacuum-polarization function. A third modification is that we incorporate the very 
recent results of Ref. M concerning the irreducible two-loop corrections of 0((a vn^ j w? w ) 2 ) . 
In the rest of the calculations, as we did in Ref. 0, we treat the u, d, and s quarks as 
massless but we now include terms of O (a m 2 / , where / = c,b,r, . . . . Although they 
are very small — they contribute to Ar only at the < 1 x 10~ 4 level — , their incorporation 
may facilitate detailed comparisons with calculations by other authors. 
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The plan of the paper is the following: in Section 2, we discuss the definitions of the 
basic Ws parameters e 2 {m z ) and sm 2 9 w (m z ) in the presence of the 0(aa s ) corrections, 
with emphasis on the decoupling of heavy particles. In Section 3, we incorporate the 
perturbative 0(aa s ) contributions in Af w and Af . We emphasize the important fact that 
the magnitude of the 0(aa s ) effects depends sensitively on the precise definition of m t . 
Our detailed calculations, as well as the other papers in the literature, employ the "on- 
shell" definition of m t . In the discussion we give a brief argument to indicate why this 
choice is useful and appropriate. In Section 4, we present a simple method to separate the 
residual, i.e. "non-electromagnetic", parts of Af, Af w , and Ar. In Section 5, we include 
the perturbative 0(aa s ) corrections in the calculation of Ar in the on-shell scheme, using 
the formulation of Ref. M. In Section 6, we discuss the incorporation of threshold effects 
in Af, Af w , and Ar. We include estimates, based on a simple "Bohr-atom" model, of the 
contribution of the IS toponium resonance to the imaginary part of the self energies. We 
find that this simple model gives values roughly similar to the calculations carried out with 
more realistic quark-ant iquark potentials. In Section 7, we use the theoretical results to 
carry out precise calculations of m w and sin 2 ^(m z ), in the manner explained before. We 
find that the QCD and other higher-order corrections induce shifts in Tn w and sin O w (tti z ^ 
comparable with the expected experimental precision. Interestingly, all of them increase 
the value of m t derived from current measurements. We also compare the Ws and on-shell 
calculation of Ar, in the formulation of Ref. ||. We find that, as was the case in the absence 
of perturbative 0(aa s ) corrections, the Ws and on-shell calculations of Ar are very close 
over a large range of m t and m H values. The Appendices discuss basic expressions for the 
perturbative (D(aa s ) corrections, the effect of top-quark decoupling in Af w and Af , and 
the very small contribution from finite fermion masses. 



2 The Parameters e 2 (m z ) and sm 2 9 w (m z ) 

In our previous treatment [I] we defined these parameters, at the one-loop level, by simply 
subtracting from the radiatively corrected cofactors the terms involving 

5 = -l I + i[ 7 -l n (4 7 r)] ) (5a) 

and setting the 't Hooft mass scale, /z, equal to m z . Because at one loop S always occurs 
in combination with — In /z, this is equivalent to subtracting only the pole terms, (n — 4)~ 1 , 
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rescaling fi according to 



\x e 



7 /2 



(5b) 



(4tt) 1 /2' 

and then setting /i' = m z . The second formulation can be conveniently generalized to 
higher-order corrections and one can define the Ws renormalization procedure as the sub- 
traction of pole terms of the form (n— 4)~ m , where m is an integer > 1, and the identification 
of the rescaled 't Hooft parameter \J with the relevant mass scale, in this case m z . As is well 
known, the factor e 7 / 2 (47r) _1//2 is appended in Eq. (|5b|) to cancel relatively large numerical 
constants that are an artifact of dimensional regularization [[34]]. 

In Ref. |IJ we applied this procedure uniformly, independently of whether the top quark 
is more or less massive than m z . On the other hand, it is desirable to treat heavy particles, 
as much as possible, as decoupled. For example, when m t > m z it is convenient to subtract 
from the amplitude terms involving \n(rn t /m z ) and to absorb them in the coupling con- 
stants. In Ref. |jj we did not follow this route for two reasons: i) sin 2 9 w (m z ) appears as a 
cofactor in several important radiatively corrected amplitudes and it is not possible to ab- 
sorb completely the ln(m t /m z ) terms occurring in all of them; ii) some important relations, 
such as Eq. (D), contain terms proportional to m 2 , which certainly do not decouple. Re- 
cently, however, Marciano and Rosner [0,0 proposed to implement the decoupling idea, 



at the one-loop level, by absorbing in sin 2 ^ 



m 5 



all \n(m/m z ) terms with m > m z occur- 



ring in a specific amplitude, namely the self energy evaluated at q 2 



m. 



-ReA lZ {m 2 z ). 



Here m is the mass of the top quark or any unknown heavy particle with m > m z . With 
this convention, a heavy top or a heavy unknown particle decouples in the limit m/m z ^> 1 
from the cofactor k(q 2 ) multiplying s 2 in most neutral-current processes; as a consequence, 
this parameter can effectively be determined from the on-resonance asymmetries without 
hindrance from unknown "heavy physics." 

We now explain how we implement the decoupling idea in the presence of the 0(aa s ) 
corrections. It is convenient to start with e 2 (m z ), which we frequently abbreviate as e 2 . We 
recall the relation between the bare charge eo and the conventional renormalized charge e: 



e 2 = e 2 l + e 2 n(/)(0) + 



7e 2 



1 



n - 4 



+ In 



m w 



1 

21 



-l 



(6) 



where 11^(0) is the usual fermionic vacuum-polarization function evaluated at q 2 = 



and the last term represents 0(efy bosonic contributions to charge renormalization that 
must be included in the SM. The latter can be gleaned, for example, from Ref. 0. As 
explained in Refs. (l|,|2|, because of the existence of mass singularities associated with the 
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light quarks, it is not possible to calculate perturbatively the contribution of the first five 
quark flavors to n^(0). Calling this contribution II® (0), the problem is circumvented by 
writing II® (0) = Ve (n® (0) - II® (m|)) IT® (m|). The first term represents the five- 
flavor contribution to the renormalized photon vacuum-polarization function at q 2 = m 2 
and can be evaluated using dispersion relations, experimental data on e + e~ — ► hadrons and 
QCD corrections []2"|, |9|, |i~0~| . The presence of a large invariant momentum, q 2 = m|, in the 
third term prevents the occurrence of fermionic mass singularities and, as a consequence, 
7&n®(m|) can be analysed perturbatively. The same is, of course, true for the leptonic 
and bosonic contributions and, because of its large mass, for the top contribution nj^(0). 
Including irreducible two-loop contributions of 0(dca s , d 2 ) to e 2 n^(0) arising from virtual 
gluon and photon interchanges, we obtain (see Appendix A) 

2d 



n$(o) 



3tt 



E 







45d~ 




+ 327T 


mi 
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8a 
9vr 
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VTLr, 



q+t 

+e 2 fle(n(5,»(0)-n< 5 > 



2]n^fl + | + |«! + + & + 3 iQl U(r.) 



a l 
7T (n -4)' 



(7a) 



where the I and q sums run over leptons and quarks, respectively, the color factor 3 is 
henceforth explicitly included, r q = m 2 / (4m 2 ,), 



/ 1 (r) = ln(4r)- ( 2 + - 



f 2 (r) = ln(4r) + 



3 

2r V r' 
KeVAr 



IX 1/2 

1 ) cosh 1 



55 



4C(3)+ 10 . 
r 12 

55 3 

4C(3) + ---ln(4r ^ _ 



5 1 

r + - + ~ 

3 r 



(7b) 



(7c) 



C(3) = 1.20206. . ., Vi(r) is a complicated function defined in Ref. [p5fl , and / = 16/3 + 
5d s /(37r) + lld/(97r). The first and second terms in Eq. ( |7a| ) are the finite parts of the 
leptonic and top contributions, while the sum over q is the finite part of the perturbative 
evaluation of e 2 Tie II® (m 2 z ). In the latter the terms proportional to r~ l are extremely small 
and we can replace the functions fi(r) and ^(t*) by their asymptotic values /i(oo) = 5/3 
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and /2(oo) = 55/12 — 4^(3) ~ —0.22491. The last term in Eq. flTa] ) represents the divergent 
part up to terms of 0{& 2 , <M S ). As is well known, up to two loops Il 77 (g 2 ) is linear in 
(n — 4) _1 and In //. We note that the cofactors of (n — 4) _1 are equal to those of ln(l///) 
at one loop and to one-half those of ln(l///) at two loops. 

In order to obtain the relation between e 2 and e 2 , one writes = e 2 /Z e in Eq. (j^), uses 
the counterterms present in Z e to cancel the (n — 4)~ x terms in Eqs. (|6]) and (|7a|), and sets 
/i' = m z in those equations. The mass scale employed in a s is discussed in greater detail at 
the end of this section. With the exception of the top quark, all the particles contributing 
to Eqs. (||[7aD are less massive than m z and their contribution is retained. To implement 
decoupling in Eq. (^]), we also subtract the finite top contribution when m t > m z , so that 
Z e in that case contains an additional finite counterterm and reads 



7T V 2 



n 



e 2 Um)_9(m t 



m 



(7d) 



where ms denotes the "finite part" after the ms renormalization has been carried out, i.e. 
the remainder after the (n — 4) poles have been subtracted and // has been set equal to 



m 



z- 



and the superscript (t) refers to the top-quark contribution. Specifically, 



n«(o) 



9tt 




(7e) 



The term proportional to (n — 4) _1 in Eq. (|73] ) cancels the divergent parts in Eqs 
The other term subtracts the finite top contribution, i.e. the second term in Eq. (|7. 
when m t > m z . Because this contribution does not exactly vanish at m t = m 



the above 



prescription leads to a small discontinuity at m t = m z . An alternative that would ensure 
continuity would be to do the matching at the point where the finite top correction vanishes, 
which is m t = 1.073 m z = 97.8 GeV. However, because the decoupling at m t > m z is easy 



to implement and is analogous to what is done in some QCD calculations [|3~5[1 , we will 
adopt it as our convention. 

As the mass range m t < 91 GeV has been excluded at the 95% confidence level |JB 
will henceforth assume that m t > m z , in which case Eqs. (|6|,[7a|,d) lead to 



we 



1 + («/7T)A 



3a) 
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where 




An 2 -Re (n$(0)-ng( 



(8b) 



and c 2 is an abbreviation for cos 2 6 W = m^/m 2 . Solving for e 2 , we have 



1 - (a/7r)A„ 



ic) 



In (a/7r)A 7 we have retained very small terms of 0(a 2 ) arising from virtual-photon inter- 
change. We have done this because of the analogy with gluon contributions and the fact 
that they contain interesting leptonic mass singularities. It should be understood, however, 
that this does not represent a complete 0(a 2 ) calculation as there are other irreducible two- 
loop contributions of this order arising from the bosonic sector and from W^, Z°, and H 
interchanges in the fermionic sector. 

Equation fl8cj ) allows us to compute e 2 in terms of a, independently of m t or un- 
known particles heavier than m z . Using m z = 91.187 GeV [fL6fl, m w = 80.22 GeV fl37[,|3g 



a s = a s {m z ) = 0.118 ||, and e 2 Tie (lljg(0) - Ii^{m 2 z )) = 0.0282 ± 0.0009 0, we find 
(«/tt)A 7 = 0.0668 ± 0.0009 or a' 1 = (4vr/e 2 ) = 127.9 ± 0.1. The expression for (a/vr)A 7 
differs from the quantity — 25e/e\jjs m Ref. |Q by the exclusion of the ]n(mt/m z ) terms, 
a more accurate description of the QCD corrections, and the updated calculation of the 
five-flavor contribution. Numerically, however, d -1 is very close to the value 127.8 ± 0.1 



reported in Ref. ffl for m t 



m 



z- 



the small difference essentially arising from the change in 



the five-flavor contribution ||10|| . 

Concerning s 2 = sm 2 9 w (m z ), we recall that in the neutral-current amplitudes this 
parameter is multiplied by the electroweak form factor k, which contains the 'jZ mixing 
term —(c/s)A lZ (q 2 )/q 2 |17l, [40|l . Here A lZ (q 2 ) is the unrenormalized transverse mixing 
amplitude as defined in Ref. expressed in terms of the Is couplings e, s, and c. In order 



to implement the decoupling, we apply the Marciano-Rosner convention | 32| , |33|| , according 
to which the ln(m f /m z ) terms in Tie A lZ {m 2 z )/m 2 z are subtracted for m t > m z . At the 
two-loop level there is also an mrindependent term, which must be subtracted, too. More 
generally, the idea is to subtract all contributions to Tie A lZ (m 2 z ) / m 2 z that involve particles 
of mass m > m z and do not decouple in the limit m — > oo. This can be implemented 
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by adding a finite counterterm in the s 2 renormalization, in analogy with Eq. (|7d|). Up to 
terms of 0(aa s ) we find, in the case of the top quark, 



1 + 



cte c 

7 + ~ 

n — 4 s 



(9a) 



and ms and the super- 



where Sq is the bare parameter, (A^z) (0) 
script (£) have the same meaning as in Eq. (|7d[) . Again, cte/(n — 4) is the divergent part of 
the counterterm and the last term is a finite contribution necessary to implement the de- 
coupling of the ln(mt/m^) and constant terms in the top contribution to Ve A lZ {m 2 z )/rn 2 z . 
Specifically, we find: 

(9b) 



K'i)'e» 



a 



where 



d 



In 



m t 



1 + 



o:. 



15cL 



7T 



■S7T 



(9c) 



In Eqs. (|9a|— c) we have neglected all terms of 0(a 2 ). 

We must still discuss the mass scale employed in a s . Following Refs. f22|, p4|-|26|, pl| , in 
the present paper corrections arising from the (t,b) isodoublet are computed with ct s (mt). 
The reason is that the dominant contributions due to this isodoublet involve mass scales 
of 0{m t ) pi|,^2[. This choice can also be justified by arguments based on effective field 
theory |]43 |. On the other hand, the perturbative contributions from the two light-quark 
isodoublets involve self energies evaluated at q 2 = rri^ or q 2 = m 2 and, for that reason, are 
calculated with a(m z ). As an example, in Eqs. (|7|,^c[), which involve top-quark contribu- 
tions, we identify a s with a s (m t ), while in the perturbative part of Eq. (jSbD , which includes 
light-quark isodoublets, we employ a s {m z ). (In the latter equation we make a slight and 
negligible change to the above isodoublet rule by also evaluating the very small bottom 
contribution with a s (m z ).) Numerically, the finite counterterms in Eqs. (|75| , |9a|) are quite 
small: for m t = 150 GeV, they are 1.0 x 10~ 3 in Eqs. flTd]-e) and 6.1 x 10~ 4 in Eqs. (|9~5|-c); 
for m t = 250 GeV, the corresponding values are 2.2 x 10~ 3 and 1.5 x 10~ 3 , respectively. 
As we will see, when s 2 is defined according to the decoupling convention explained above, 
the finite counterterm in Eq. (£Ja|) introduces small compensatory shifts in Af w and Af. 
Similarly, Z e , defined in Eq. (|7d|) , will introduce small compensatory changes in radiative 
corrections whenever e 2 is employed as zeroth-order parameter. 
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3 Perturbative 0(aa s ) Corrections to Ar w , Ar, and Ar 



When the decoupling of the top quark is implemented according to the discussion of Sec- 
tion 2, the expression for Af w (cf. Eq.(Q)) becomes (see Appendix B) 



An 




{j^w) A ww 



(0)' 



mt 



16tt 2 s 2 



6 + 



In c 2 



2~2 S 
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where A ww (q 2 ) is the unrenormalized WW transverse self energy with the coupling e 2 /s 2 
factored out ffl, d is defined in Eq. (0), and the Ms symbol has the same meaning as 
in Section 2. The quantities Ar w , Af, and Ar are gauge-invariant but some of their 
partial bosonic contributions are not. In Eq. ( |10a| ) and henceforth all explicit expressions 
and partial contributions are given in the 't Hooft-Feynman gauge. In particular, the 
expression involving curly brackets represents vertex- and box-diagram corrections to fi 
decay, evaluated in that gauge. 

Except for the last term, Eq. ( |10a| ) is the same as Eq. (8b) in Ref . |1| . As s 2 enters Eq. (Q) 
as a zeroth-order parameter and is defined in the present paper according to the decoupling 
convention explained in Section 2, there is now an additional contribution — (a/ir)d arising 
from the finite contribution in Eqs. (|9a]-c). 



We now discuss the evaluation of Eq. ( |10a| ) . The quantity 
from Eq. ([7a]). We have 
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where A 7 is defined in Eq. ( [Sb]) and the other terms on the r.h.s. represent the top con- 
tribution. Inserting Eq. (|10b| ) into Eq. ( |10a| ) and neglecting small terms of 0(a 2 ), we 
obtain 
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where 



d 

l-8s 2 /3' 



(lOd) 



The bosonic contribution (e 2 / s 2 m1 v ) [Aww(n^w)— Aww(fy]ws is given in Eq. (A. 6) of Ref. 
In order to study the fermionic contributions, denoted by a superscript (/), we define 



B (f) 



-Re 
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and write 



where Bq and -Bqc D stand for the 0(a) and 0(aa s ) corrections. If we neglect 
terms proportional to ami/m^ (q = d,s,b), the mixing angles in the quark 
irrelevant 



BW = B { n f) + B 
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very small 
sector are 
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and to zeroth order in a s we have (cf. Eq. (A. 5) of Ref. M 

I) 2 
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where we have included the lepton and quark contributions and u t = m 2 /m 2 v . As a 
refinement, in Appendix C we discuss the contributions of (D(m 2 /m^), where m stands 
for a lepton or quark mass other than m t . As these terms are very small, we may neglect 
in their evaluation the squares of the mixing angles. In that case one obtains a sum of 
isodoublet contributions, which, for arbitrary masses, is given in Appendix C. One finds 
that the corrections of 0(m 2 /m^) to Bq^ are indeed very small, of O(l0~ 5 ). 
The contribution of 0(aa s ) in Eq. ( |llbp is given by (see Appendix A): 



B, 



(/) 

QCD 
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where F±(x) is defined in Ref. [^5| . In Eq. ( |lld|) we have neglected all quark masses other 
than m t . The first term in Eq. ( |lld|) arises from the (u, d) and (c, s) isodoublets, while the 
second is associated with the (t, b) doublet. In Eqs. ( |10c| -d, pTd| ) the mass scale of a s has 
been chosen according to the prescription explained at the end of Section 2. 

The asymptotic behavior of Af w for large m t can be obtained from Eqs. (]9c , pL 0c| — d. , |1 1 cj — 
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d) and is given by 



a * a / 1 , a s{ m t) \ , m 

Ar w ~ —2 1 + In 

07TS \ 7T / m z 



(m t > m z ). 



;i2) 



This exhibits a smaller coefficient than the corresponding expression in Ref. ||12|| , a feature 
that is due to a partial cancellation with the finite counterterm — (a/7r)d in Eq. (|10a|) . The 
asymptotic behavior for m H ^> m z is the same as in Ref. ||12|| , namely 
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We now turn our attention to Af. When the top decoupling is implemented according 
to Section 2, the expression for Af becomes (see Appendix B) 



Af = Afar — 
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where o? and d are defined in Eqs. (|9cl , |10d| ) and A zz (q 2 ) is the unrenormalized ZZ transverse 
self energy with the coupling e 2 /s 2 factored out [TJ. Except for the terms involving d and 
d, which arise from the finite counterterms in Eqs. ( ffd] , |9aD associated with the decoupling 
of the top quark, Eq. (14) has the same form as Eq. (15b) of Ref. [[jj]. It is understood, 
however, that Af w , e 2 , and s 2 in Eq. (|TJj) are defined according to the prescriptions of the 
present paper, namely Eqs. ( |10c| , |3c| ,|2"|). In Appendix B we show how Eqs. ( |10a| -c,|14"D can 
be derived from the results of Ref. [|TJ by neglecting very small contributions of O (a 2 ) 
without logarithmic or m 2 /m 2 enhancements, as well as terms of 0(a 3 ). As pointed 
out in Ref. |J, if one neglects also subleading corrections of 0((a/irs 2 )xt), with x t being 
defined in Eq. ( |17b|) , one can replace e 2 1 — Ar w — (a/n)d — > e 2 in the second term of 
Eq. flT4"p. As in Ref. ||, we have retained such subleading terms in Eq. flTj] ) because the 
resulting expression describes very accurately the resummation of one-loop effects, which 
is particularly simple in the ms framework. 

We now turn our attention to the evaluation of the self energies in the second term of 



Eq. fli~4]). The bosonic contribution (e 2 / s 2 m 2 z ) Tie A\J w (m 2 w ) / c 2 — A zz (m 2 ) 
Eq. (A. 9) of Ref. [p]]. To study the fermionic contributions, we write 
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and expand 



where Cq and Cqq D are the O(a) and 0(aa s ) corrections, respectively. If the small 
terms proportional to a m^/m^, (q = d, s, b) are neglected, the mixing angles are once more 
irrelevant and, to zeroth order in a s , we have (cf. Eq. (A. 8) of Ref. [|TJ) 
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where /j, t = m 2 /m 2 z , D t = Afi t — 1, and A(D) = D 1 / 2 tan 1 (D 1 I 2 ) for £> > 0. In analogy 

if) ( f) 

with our discussion of Bq , in Appendix C we give the expression for C for arbitrary 

fermion masses in the approximation of neglecting the squares of mixing angles. We find 

that the corrections of 0(m 2 /m 2 ^-) to Cq are also of 0(1O~ 5 ). 

The contribution of 0(aa s ) in Eq. ( |15bp is given by 
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where the functions Fi(x), Vi(r), and Ai(r) are defined in Ref. [p5|1 . 

The leading asymptotic behavior of Cq for large m t is (3a/167rs 2 c 2 )(m 2 /m|), which 
arises from the last term of Eq. ( |15c| ). In the 0(aa s ) corrections the leading contribu- 
tion is contained in the combination (aa s (m t ) /47r 2 s 2 c 2 ) |// t Ai(l/(4// t )) — 4c 2 co , t F 1 (l/c<j t )], 
which asymptotically becomes — (aa s (m t )/S7T 2 s 2 c 2 ) (tt 2 /3 + 1) (mf/m 2 ). Combining these 
contributions, inserting the result in Eq. ([14]) and neglecting there the subleading 0(a 2 ) 
contributions, we have 
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while the leading asymptotic behavior for large m H is the same as in Ref. namely 



Oi ( 5 3 \ TYL 

Af ~ - - -c 2 In— (m H >m z ). (16b) 

27^^ \6 4 / m z 

The term involving d s in the last factor of Eq. ( |16a|) represents the most important 0(aa s ) 
correction. Indeed, the contribution of Eq. (|16a|) depends very sensitively on m t and the 
coefficient of a s /ir in the last factor, namely —2.86, is quite large. The presence of this a s 
correction induces an increase in the value of m t read from experiments of approximately 

Am t a s {m t ) fn 2 \ 

— - — + 0.455 a a (mt), (16c) 

m t 3n \ 3 J 

which amounts to Am t « (4.9, 7.5, 9.6, 11.7) GeV for m t = m z , 150, 200, 250 GeV. It 
should be stressed that, as is obvious from the structure of Eq. (|16a|) , these results depend 
sensitively on the precise definition of m t . The quantity that appears in Eq. (|16a|) and the 
various expressions of this paper is the zero of the real part of the inverse propagator. In 
the literature it is variously referred to as the "physical," "on-shell," or "dressed" mass. In 
the approximation of neglecting the s dependence of the imaginary part of the top-quark 



self energy, it coincides with the real part of the complex pole position |[T3Hl5[1. It is also 
the mass that occurs in the Balmer formula for the toponium levels in the non-relativistic 
bound-state picture and the parameter that governs the start of the ti cut in perturbation 



theory |44j]. All the recent calculations of 0(aa s ) and tt production p2"|-j3T|l employ this def- 



inition or slight modifications thereof. It is worthwhile to notice that the 0(a s ) corrections 
become much smaller if one employs other definitions of mass |45j . For example, m t and the 



41 . Insert- 



running mass, m t (m 4 ), are related by m t = m 4 (m t ) 1 + 4d s /(37r) + Oya 
ing this into Eq. ( |16a| ), we get a contribution involving (rh t (rh t )) 2 [1 — 2(d s /7r)(7r 2 /9 — 1)] 
and we see that the coefficient of a s /n has changed from —2.86 to —0.19. Similarly, if 
one expresses the corrections in terms of the Georgi-Politzer mass, M(—m 2 ) |^6|], which 



is gauge-dependent and usually evaluated in the Landau gauge, the coefficient of a s /ii 
becomes even smaller, namely +0.09. Because the perturbative evaluation of the (t, b)- 
isodoublet loops involves high mass scales, of 0(mt) or 0(m z ), both the on-shell and 
rht{rht) definitions are in principle suitable, although the former is the natural choice in the 
DR approach p3Vf2B],pil|,|3l|. The relevant question, of course, is what mass parametriza- 
tion is more adequate to describe the physical issues at hand, namely the production and 
detection of the top quark. In this connection it also appears that the on-shell mass is the 
most appropriate parameter because, in the propagation of t and i between the "produc- 
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tion" and "decay" vertices, configurations near the "mass-shell" will be greatly enhanced 
kinematically (resonance effect). Another consideration, of a more practical nature, is that 
the writ parametrization of the radiative corrections should be consistent with the one em- 
ployed in the calculation of tt production |?7j-29l and, as mentioned before, this is again 
the "pole" or "on-shell" definition. 

Returning to the evaluation of Af , we must still consider the irreducible contributions of 
0{a 2 {m 2 /m 2 v ) 2 ). As mentioned before, these can be gleaned from the two-loop irreducible 
corrections to the p parameter. The very recent work of Ref. ]7| leads to a significant 
change in the magnitude of these corrections. Indeed, these authors find for the leading 
high-mj contributions to the p parameter an expansion of the form 

r 1 -*^*©!) 31 -** (17a) 

where 

x t = (17b) 



is the one-loop term [47J and R is a negative function of m H /m t . When m H /m t = 0, R 
equals 19— 2n 2 rs —0.7392, the result of Ref. |3|, but as m H /m t increases, R rapidly becomes 
more negative, reaching a minimum of ~ —11.8 for m H /m t ~ 5.8. In current discussions 
m H /m t ranges from « 0.24 (corresponding to m H ~ 60 GeV and m t ~ 250 GeV) to ~ 11 
(corresponding to m H « 1 TeV and m t « 91 GeV). Noting that i?(0.24) w —3 and 
R(ll) ~ —10, it is clear that, although Refs. and agree in the limit m H /m t = 0, 
for realistic values of this ratio the results of Ref. tell us that these corrections are 
considerably larger in magnitude. 

Calling Af^ the "one-loop" expression for Af given in Eq. fll4]), we include the two-loop 
irreducible contributions by writing 

Af« x2 



Af = Af« --R(^-) x\ ^ — (18a) 

The rationale is the following. In the Ws scheme the p parameter is naturally identified 
with p = c 2 /c 2 and, from Eqs. (0J|), we see that (cf. Eqs. (17a,b) of Ref. [I|]) 

1 — Af 

Neglecting very small terms of order (Rxf/3) 2 , one indeed verifies that when the sec- 
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ond term in Eq. ( |18a| ) is inserted in Eq. ( |18b| ), it leads to an additional contribution of 
—R{m H /m t )x 2 /?> to 1/p, in conformity with Eq. ( |17a| ). 

We have given all the elements necessary to evaluate the basic radiative corrections Ar w 
(cf. Eq. ( |10c|) ) and Af (cf. Eqs. (|l4 i|18af )) including 0(aa s ) corrections. In conjunction with 



AA 2 



nl/2' 



m?Jl — Ar) 



19a) 



which follows from Eq. (j^), Af can be employed to calculate s 2 = sm 2 9 w (m z ) in terms 
of the accurately known quantities G^, m z , and a, as a function of m t and m H . The 
parameter s 2 = 1 — m^/m 2 can be computed from (cf. Eq. (19) of Ref. M) 



1 - T7T 



c 2 Af w — Af 



1 - Af, 



(19b) 



and Ar (cf. Eq. (22) of Ref. |]) from 

;,2 



Ar = Af 
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w — ~ 2 



Af w — Af 



s 2 1 - (c 2 /s 2 ) (Af w - Af ) / (1 - Arv) ' 
Alternatively, writing p = 1/(1 — Ap) we have (cf. Eqs. (17a,b,20) of Ref. M) 

Ar w - Af 



Ap 
1 - Ar 



Af 



r Ap )(1-Ar w ) 



(19c) 



(19d) 
(19e) 



Again, the m w -m z interdependence can be expressed in two equivalent forms (cf. Eqs. (24,25) 
of Ref. mi), 



,2 



2< 1 + 
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m 2 (l — Ar) 



ml 2 



4A 2 
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m 2 z p(l - Af w ) 



(19f) 
(19g) 



The mass can be evaluated from ( |19b|) , or (|19f|) , or (|19g|) . Equations ( |19a| -g) have the 
same structure as in Ref. [jXJ] because they follow from the same basic relations, namely 
Eqs. (|l|-0). In this paper, however, the explicit evaluation of Af w and Af via Eqs. 



|18a|) is somewhat different because we have included the 0(aa s ) corrections, updated 
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the contributions of e 2 Tie (n^(0) — n^(m|)J and the two-loop irreducible parts, and 
implemented the decoupling of the top quark. By inference the same holds true for Ar and 
Ap when they are evaluated from Af w and Af via Eqs. ( |19c|) and (|19d| ), respectively. It 
should be observed moreover that, except for very small effects of 0(a 2 ), i.e. of the same 
order as those we have neglected, the decoupling of the top quark should not affect physical 
observables such as the radiative correction Ar and the predicted value of m w . 

As shown in Appendix C, the corrections of 0(m 2 j/m^j to Ar are also very small. 
However, they are enhanced relative to those of Ar w and Af, and for m z < m t < 250 GeV 
they vary from ps — 7 x 1CT 5 to « 8 x 1CT 5 . 



4 Residual Parts of Af w , Ar, and Ar 

It is a simple matter to derive expressions relating m w and m z to G^, s 2 , c 2 , and e 2 (rather 
than e 2 ). To see this, we write Eq. (H) in the form 

V2G,m 2 w 1 - (a/ir) A 7 - (Ar w - {a fir) A 7 ) ' 

Factoring out (1 — (a/7r)A 7 ) and recalling Eq. (pcj), we have 

~o 7rci 1 . , „ 

V2G^m 2 w 1 - (Ar w ) res 



where 



a, \e 2 



7 ; 2 
7i / e z 



(Arv) rcs =[Af w - -A 7 -. (20c) 



The correction (Ar w ) Tes represents the "residual part" of Af w after we have subtracted the 
large contribution (a/7r)A 7 , evaluated with the coupling e 2 rather than e 2 . As (Af w ) TCS <^ 
Ar w , we see that the dominant part of Af w can be absorbed by employing a rather than 
a as zeroth-order coupling. 

Starting with Eq. @, the analogous argument leads to 

s 2 c 2 = r Ka - , (21a) 

V^G^ml l-(Af) res ' 1 ; 

where 

(Af) rcs = (Af - ^A 7 ) J. (21b) 
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Again (Af) res is the "residual part" of Af . In connection with the inclusion of the two- 
loop irreducible contributions of 0(a 2 (m^/m^ v ) 2 ) (cf. Eq. ( |18a| )), it is easy to see that the 
one- loop" (A f^J can be obtained by replacing Af w — > (Af w ) Tes , (a/n)d — > (d/n)d, 
and {a/n)d — > (a/n)d on the r.h.s. of Eq. (0), and the final (Af) ros follows by substituting 
AfC 1 ) — > ^Af^J and Af w — > (Ar w ) Tes on the r.h.s. of Eq. (|18a| ). Of course, once Af is 
known, one can directly use Eq. (|21b|) for the numerical evaluation of (Af ) res . 
The corresponding expression from Eq. ([[]) is 

9 7TQ! 1 \ 

s 2 = __ - ——, (22a) 



V2G,m 2 v 1 - (Ar; 



res 



where 



a „ \ e 2 



7T '/ e 2 



(Ar) rcs = Ar--A 7 - (22b) 



is the "residual part" of Ar. One readily finds that this quantity can be obtained by simply 
substituting Af w — ► (Af w ) res and Af — > (Af) res everywhere on the r.h.s. of Eq. (|19cj) . Of 
course, Eq. ( |22b| ) can be directly used for numerical evaluations. 

As illustrations, for m H = 250 GeV and m t = 150 GeV we have: Af w = 7.02 x 10 -2 , 
Af = 6.34 x 10" 2 , and Ar = 4.74 x 10~ 2 , while (Af w ) rcs = 3.6 x 10~ 3 , (Af) rcs = -3.6 x 
10~ 3 , and (Ar) res = —2.08 x 10~ 2 . The corresponding values for m H = 250 GeV and 
m t = 200 GeV are Af w = 7.08 x 10~ 2 , Af = 5.88 x 10~ 2 , and Ar = 2.90 x 10" 2 , while 
(Ar w ) res = 4.3 x 10- 3 , (Af) rcs = -8.5 x 10~ 3 , and (Ar) rcs = -4.05 x 10~ 2 . Unlike {Af w ) ves 
or (Af) res , (Ar) res is quite large for m t ~ 200 GeV. 



A correction similar to (Ar) res has been recently employed in Ref. |20]. The two quanti- 
ties are, however, not identical because a, defined in the Ws scheme, differs somewhat from 
the effective parameter a(m z ) = (128. 8)" 1 used in that work. This illustrates the rather 
obvious but important fact that running couplings are scheme-dependent. 



5 Ar in the On-Shell Scheme 

In this section we discuss the incorporation of the perturbative 0(aa s ) corrections to 
Ar [0,0] in the on-shell scheme of renormalization 0. We follow the formulation proposed 
recently in Ref. ||, based on the expression 

c 2 

Ar = Aa - —Ap (1 - Aa) + (Ar) rem , (23a) 
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where Aa = 0.0597 ± 0.0009 represents the contribution of the charged leptons and the 
first five quark flavors to the photon vacuum-polarization function evaluated at q 2 = m|, 
i.e. e 2 73e (IT 77 (0) - IT 77 (m|)), 

Ap = x t (l + R (^) f) , (23b) 

x t is defined in Eq. (|17b| ), and R{m H /m t ) [|7| is the function discussed after Eq. ( |17a| ). The 
result for Aa quoted above includes the recent calculation of the first five quark flavor 
contribution [K| and, for this reason, it slightly differs from the central value of 0.0602 
employed in Ref. ||. The second term in Eq. ( |23a|) involves the leading m t - dependent 



correction Ap to 1 — 1/p (cf. Eq. (|17a|) ) and we see that, in the case of Ar, it is enhanced 
by a factor c 2 /s 2 . Its importance for large m t in the m w -m z interdependence was pointed 
out in 1980, in the work of W. J. Marciano and one of us (A.S.) ||41|| . Since that time, this 
potential effect has been discussed by several authors H,|5|,p2T^6|, ^8l[49| . 



The separation into leading contributions (the first two terms in Eq. ( |23a| )) and a 
'remainder" (Ar) rem is the same as was proposed in Refs. ||,|5|,[26|], except that we have 



included the recent results of Ref. on the two-loop contribution to Ap. It is important 
to note that (Ar) rem differs from the quantity (Ar) res introduced at the end of Section 3. 
Whereas in the latter we subtract the large logarithmic corrections, in (Ar) rem we also 
exclude the leading m r dependent contributions. The formulation of Ref. |§ provides also 
a very specific prescription to calculate (Ar) rem , namely 

i * \ . t i\ . c 2 ^ c 2 V2G u m 2 (1 — Aa)s 2 
Ar rcm = Ar« -Aa + -X + -{x t - X)- ^-^ J —, 24a 

where Ar^ is the familiar one-loop expression of Ref. 0, expressed in terms of a and a/s 2 
as expansion parameters, 

3a m 2 (c>av.\ 
Xt 16tcs 2 m 2 w ' 



e 2 

X = —Re 

s 2 



A ww (m 2 v ) A zz (m 2 z ) 



m 2 w m\ 



(24c) 

fin 



and the subscript fin means "finite part", i.e. that the pole terms have been subtracted 
and \j! has been set equal to m z . We also note that X is a gauge-invariant quantity. In 
Eq. (|24cj) we follow the notation of Ref. [Q, which differs from that of Refs. |2, 41] in that 



an explicit coupling e 2 /s 2 has been factored out in the A ww and A zz self energies. The 
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term — Aa in Eq. ( |24a ) subtracts from Ar^ the large logarithmic corrections which are 
included as part of the leading contributions in Eq. ( |23a| ). The quantity (c 2 /s 2 )X subtracts 
another important part of Ar^\ which is then treated more accurately in the second term 
of Eq. ( 23a ) and the third term of Eq. ( p4a| ) , according to the following prescription. 
Decomposing — (c 2 /s 2 )X = — (c 2 /s 2 )x t + (c 2 /s 2 )(x t — X), the dominant part, — (c 2 /s 2 )x t , 
is included in the second term of Eq. (|23a|) with the effective coupling constant changed 
according to a/s 2 — > \[2G il w? w {\ — Aa)/7r. The non-dominant part, (c 2 /s 2 )(x t — X), is 
treated with the same coupling modification but it is retained as part of (Ar) rem (third term 
of Eq. (|24a|) ) . The rationale for this treatment of (c 2 /s 2 )(x t — X) was explained in Refs. 
[P|,|50| and reflects the fact that a careful analysis of the resummation of one-loop effects leads 
to an expression of the form of Eq. ( |23a| ) in which X, rather than its dominant part Xt, is 
multiplied by ^/2G^ t s 2 m 2 v (l—Aa)/(Tca) = (1— Aa)/(1— Ar). This is natural because, since 
Ar is the radiative correction in the relation between m w , m z , G^, and a, it should involve 
quantities evaluated at q 2 = m 2 and q 2 = such as X, rather than amplitudes evaluated 
at q 2 = such as Ap. We also point out that the neglect of this effect in Eq. (|24a|) , namely 
the replacement v / 2G M m^s 2 (l — Aa)/(na) — > 1 in the last term, would induce a change 
(c 2 /s 2 )(x t — X)(c 2 /s 2 )Ap = 0((a/irs 2 )(c 2 / s 2 ) 2 Ap)] although formally subleading, this 
is enhanced by two powers of c 2 / s 2 and is, therefore, significantly larger than the expected 
theoretical error. It was already pointed out in Ref. || that Eqs. ( |23a| , |24a| ) include correctly 
not only the leading terms of 0{a 2 {m 2 /m 2 v ) 2 ) and 0{a 2 ln 2 (m z /m/)) , where rrif is a 
generic fermion mass, but also the subleading contributions of 0(a 2 ln(m z /rrif)). 

We now turn our attention to the incorporation of the perturbative 0(aa s ) contri- 
butions. We first consider Eq. ([24a|) . The amplitudes modified by the QCD correc- 



tions are X"> in the last term and the self-energy contributions (e 2 /s 2 m 2 v )Re 



—Aww(0) and e 2 IIj$(0) contained in Ar^\ where the superscript (/) denotes again 



fin 



~ry V y Jfin 



fermionic contributions. The first two are obtained from C^> (cf. Eqs. ( |15a| -d)) and 
(cf. Eqs. ( |lla| -d)), respectively, by simply changing everywhere s 2 — > s 2 and c 2 —> c 2 . In 
particular, their QCD corrections are derived from Eqs. ( |15d|) and ( |lld|) , respectively. The 
correction e 2 11^(0) can be read from Eq. ( |7a|) by removing the pole terms, setting 
pi = m z , and substituting a — > a. It is not necessary to consider (c 2 / ' s 2 )X^ in the third 
term of Eq. (|24a|) because, as explained before, it cancels an identical contribution in Ar^\ 
The value of Aa is not modified, as QCD corrections have already been taken into account 
in its evaluation. 



In the above discussion, the quantities Xt in Eq. (|24a|) and Ap in Eq. ( |23a| ) have not 
been altered, so that, except for Aa, all the QCD corrections are contained in (Ar) rem . 
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We may wish, however, to incorporate the leading QCD corrections in the second term of 
Eq. ( |23a| ). To achieve this, we subtract them from Eq. Q24a| ) by replacing 



x t 



3a m 2 



16ns 2 ot?„ 



2a s {m t ) ( tt 2 
3tt 1 3 



+ 1 



(25a) 



in the last term of that equation and, at the same time, we substitute 



Ap -> Ap' 



x t 



2a s (m t ) fir 2 \ „ [m H \ x t 

-Sr 1 t + 1 ) +R {^)i 



(25b) 



in Eqs. ( |23a| ,b). The overall evaluation of Ar is, of course, the same whether we employ 
Ap in Eq. ( |23a| ) and x t in Eq. ( |24a|) or the modified quantities, Ap' and x' t . In the second 
formulation, however, the first two terms of Eq. (|23a|) , with Ap — > Ap', describe more 
accurately the leading m t -dependent corrections. 

The contributions of 0(a 2 a S} acaf) are unknown at the present time and for this reason 
we have not made any attempt to include them. However, the structure of Eq. fl25b|) gives 
a hint about what their magnitude might be. Suppose, for example, that the leading QCD 
effects are always very small when the electroweak corrections are expressed in terms of the 
running mass, m t (m t ), as it happens with the 0(aa s ) corrections. In that hypothetical 
case, the discussion after Eq. (|i6c|) indicates that the modified Ap parameter would be 



obtained approximately by appending a factor [1 — (2d s (m t )/37r) (vr 2 /3 + 1)] to each x t in 
Eq. (|23bD . For m H = 600 GeV and m t = 200 GeV, the difference with Eq. (gF| ) would lead 
to an aditional contribution to Ar of —3.4 x 10~ 4 . This is of the same order of magnitude 
as the subleading terms of (9((a/7rs 2 )(c 2 /s 2 )x 4 ), discussed in Section 1. The surprisingly 
large size of these possible corrections of C((c 2 / s 2 )a s x 2 ) is due to the m\ dependence and 
the considerable magnitude of the function R. This observation illustrates the fact that, 
for large m t values, the theoretical error due to unknown higher-order corrections may arise 
from very different sources. 



6 Threshold Effects 

The fermionic contribution to the vacuum polarization of the intermediate vector bosons 
can be expressed in terms of the amplitudes 



Uli A ( q , mi ,m 2 ) = -iJ d 4 xe^(0 \T*[J V Ax)jy\m 0), 



(26a) 
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where T* is the covariant time-ordered product and = ^17^2 and J A = ^lj^js^ are 
the vector and axial-vector currents, constructed with the spinors fields ipi and t/> 2 endowed 
with masses mi and 1712, respectively. Thus, except for vector currents with mi = rri2, the 
conservation of the currents is explicitly broken by mass terms. In turn, the tensors 11^, 
have the well-known structure 

nj^(g, mi, m 2 ) = IL v ' A (s, mi, m 2 ) + X v ' A (s, m l5 m 2 ) q^q v , (26b) 

where, throughout this section, s = q 2 . 

Threshold effects involving the ti, tb, and bb channels can be expressed as contributions 
to the imaginary parts of the functions U V ' A and X V,A p4| , p7| -|3T|l . A number of papers [p3j,|24 



have made use of DRs to express the physically important amplitudes Ii v,A in terms of their 
imaginary parts. In Refs. PD|,PTJ two of us (B.K. and A.S.) proposed to use DRs directly 



constructed from the Ward identities. We recall the basic strategy: contracting both sides 
of Eq. ( p6a| ) with g M , one derives the relation 



Tl v,A (s, mi, m 2 ) = — sX v ' A (s, mi, m 2 ) + A v,A (s, mi, m 2 ), (27a) 
where A v,A (s, mi, m 2 ) is defined by 

d 4 xe iqx (0 \T[d»jV' A (x)jV> A \0)]\ 0) = A^(s,mi,m 2 )g,. (27b) 
The idea then is to write DRs for X V ' A and A V ' A and to obtain U V ' A by means of Eq. ( |27a| ). 



This approach has been employed in Refs. [30,31] to discuss both the perturbative and the 



threshold contributions in the on-shell scheme of renormalization. The aim of this section 
is to extend the discussion, so that the threshold effects can also be included in the ms 
calculations. 

We recall that in our analysis the threshold contributions to the imaginary parts are 
non-zero over a small, finite range, i.e. they have compact support. As a consequence, the 
corresponding unsubtracted DR integrals for X V,A and A V ' A are convergent and moreover 
vanish as \s\ — > oo. In using the DR approach, we self-consistently assume that the thresh- 
old contributions to X V ' A and A V ' A tend to zero as \s\ —> oo, so that these quantities satisfy 
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unsubtracted DRs. Thus, 



X v ' (s, mi,m 2 ) 
A v ' A (s,mi,m 2 ) 



1 f , . Xm X v ' A (s', mi, m 2 ) 

— / ds ■ . 



7T 



ie 



1 /" , r Xm A v ' A (s' ,mi,m 2 ) 



7T 



s' — s — ie 



(28a) 
(28b) 



Inserting these expressions into Eq. ( |27a| ) and using the same equation to relate the imag- 
inary parts, one obtains the two equivalent representations [31] 



n 



s,m 1 ,m 2 ) 



- I ds' 

7T 



Xm II ' (s 1 , mi, 777.2) 



s' — s — it 



+ Xm X v,A (s', mi, m 2 ) 



(29a) 



Yl v ' A (s,mi,m 2 ) 



s r ds' XmIl v ' A (s',mi,m 2 ) 1 r ds' . v A , , , . , , 

- / — , 7 + - / — ImA^ s>i,m 2 . 29b 

71 J s' s — s — l€ 71 J s 



In Eqs. (p8a| - p9bD and henceforth it is understood that \ V ' A , A V ' A , and U V ' A represent 
the excess threshold contributions relative to the perturbative 0(aa s ) corrections. As 
explained in Ref. |j31| , Eq. ( |29b|) can be directly derived from the following assumptions: 
i) U V ' A satisfies a once-subtracted DR; ii) the subtraction constant is determined from the 
Ward identity (g7p, so that Tl v ' A (0) = A V > A (0); Hi) A v ' A (s) satisfies an unsubtracted DR, 
so that A V,A (0) can be calculated from the second integral in Eq. ( [29b|) . The fact that 
Y[V,a mus ^ sa tisfy a subtracted DR can be clearly seen by considering the particular case 

—U v (s',m,m)/s' and 



of vector currents with equal masses. In that case X v (s' 
A v (s',m, m) = 0, so that Eqs. ( |29a| ,b) reduce to 



m, m) 



r~\z~ f \ s [ ds' XmU v (s' ,m,m) 

11 (s, m, m) 



(30a) 



which vanishes at s = 0, in conformity with the Ward identity ( |27a| ). If, instead, H v {s, m, m) 
were to satisfy an unsubtracted DR, the condition U v (0,m,m) = would imply 
/ ds' XmU v (s' , m, m) / s' = 0, which is manifestly false, as XmU v (s' , m, m) > 0. In sum- 
mary, Eqs. (|29a|,b) are the simplest possible DRs consistent with the Ward identity 



As mentioned before, the latter is a crucial requirement. We also see from Eq. ( |29a|) that 
the threshold effects modify the asymptotic behavior of the full U v ' A (s) as \s\ — > 00 by con- 
stants, i.e. by subleading terms. (We recall that the perturbative contributions to U v ' A (s) 
grow as s, modulo logarithms.) 

Threshold effects associated with the tb channels are greatly suppressed because they 
are proportional to the squared reduced mass of the quarks [E3J] and can be neglected. 
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Those involving the bb channel, i.e. "bottomium resonances," give significant contributions 
only to e 2 Il 77 (0) and are already included in the evaluation of e 2/ Re (ll 7 ^(0) — 11^ (m|)). 
Thus, only the case m\ = m^ = m t is significant. For vector currents with equal masses, 
the relevant DR is given in Eq. fl30a|) . For axial- vector currents, we employ 



II (s,m,m) 



7T 



ds' 



Xm H A (s', m, m) 



+ Xm X A (s', m, m) 



16 



(30b) 



in accordance with Eq. fl29a| ). The amplitude XmU A (s' ,m t ,m t ) receives contributions from 
J p = 1 + states, i.e. 1 = 1, which are suppressed near threshold by centrifugal barrier effects. 
On the other hand, Xm X (s', m t , mt) also receives significant contributions from 0~ states, 
i.e. I = 0. A detailed discussion of XmH v (s',m t ,m t ) and XmX A (s',m t ,m t ) in both the 
24 1 and Green-function [28] approaches is given in Ref. f3l| . In both cases one 



resonance 



finds [01 



Xm X A (s', mt, m t ) ~ Xm X v (s', m t , m t ) 



Xm U v (s f , m t , m t ) 



(30c) 



The second equality is, of course, an exact consequence of the Ward identity ( |27a| ). 

We now discuss the specific threshold contributions to the basic radiative corrections 
studied in the paper. The corrections A 7 and Ar w contain e 2 n 7 -y(0) (cf. Eqs. (|10a| ,b)). 
The top contribution to this amplitude is obtained from Eq. ( |A.7[ ), 



e 2 nW(0) = e 2 -(n y )'(0,m 4 ,m t ), 

where (JJ v Y(0, m t , m t ) = (d/ds)U v (s, m t , m t )\ s =o- Recalling Eq. ( |30a|) , we have 

ds' 



(31a) 



1 f ds 

Cn v y(Q,mt,m t ) = - — I XmU v (s',mt,m t ) 

7T J S 



(31b) 



There are no additional significant threshold effects in Af w , because those involving 
-Re {A 

ww( m 1y) ~ A w (0)) / m w { C ^- -^l- ( |f O 3 !) ) are suppressed by reduced-mass effects. In 
the case of Af (cf. Eq. (JH|)) the same holds true for the term involving Tie A ww (mw) / (c 2 m 2 z ) ; 
there are, however, significant threshold contributions to Tie A zz {m 2 z )/m 2 z . According to 
Eq. ([A~T3D , 



Tee A$l (rn 

2 

mi 



1 

16c 2 m 2 



Tie 



8 \ 2 

1 - -s 2 ) U v (m 2 z ,m t ,m t ) + U A (m 2 z ,m t ,m t ) 



(32a) 
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Using Eqs. ( |30a| -c), we have 



Tie II y (m|, m t , m t ) 1 ^ /" rfs' 2m n v (s', m t , m t ) 



7?e II (m z , m t , m t ) 
ml 



- - v / — r ' 2 > ( 32b ) 

ir J s s' — m z 

/ ds' Xm X A (s' ,m t ,mt) 
71% J 

2 / — Trails', m t ,m t ), (32c) 



where P denotes the principal value of the integral. We note that, as explained earlier, the 
threshold behavior of Xm II (V, mt, mt) is suppressed by centrifugal barrier effects and its 
contribution has been neglected in Eq. (|32c|) . As the support for the threshold contributions 
to Xm Il y (V, m t , m t ) is located in the neighborhood of s' ~ 4m 2 , the dominant effect for 
ml m| is given by Eq. ( |32cj ) , with Eqs. ( f31b|J32b| ) being relatively suppressed by a 
factor m\j (Amf). Furthermore, the term involving TieU v (m z ,m t ,m t ) in Eq. ( |32a| ) has a 
small cofactor (1 — 8s 2 /3) 2 ~ 0.14. In the range m z < m t < 250 GeV, Eq. ( |32c| ) increases 
more rapidly than linearly with m t , while the contributions of Eqs. ( |31b| , |32bD to Af and Ar 
remain very small, at most a few times 10~ 5 . As discussed in Ref. the sign of Eq. ( |32cj ) 



can be understood with an argument reminiscent of the one employed in technicolor theories 
to explain the generation of the vector-boson masses: 0~ states contribute to the functions 
X A (s) and A A (s) and, via the Ward identity ( |27a| ), to II (s). In technicolor theories this 
gives rise to m| ~ A zz (0). In our case, when Eq. ( |32c| ) is inserted in Eq. ( |32a| ), it leads to 
a positive contribution to the Z° mass shift, 5m 2 z = Tie A zz (m 2 z ). 



Inserting into Eqs. (|31b| , |32b| ,c) the detailed evaluations of Xm U v (s f , m t , m t ) and 
Xm X A (s' ,m t ,m t ) |3T]] derived from the analyses of Refs. pi,|2^|, we obtain the additional 



contributions to Ar w and Af associated with the ti threshold. Because the threshold 
contributions in e 2 II^(0) and Tie A lZ {m 2 z ) /m 2 z are very small, we have not subtracted 
them from A 7 or the latter amplitude. Consequently, they do not affect the counterterms 
of Eqs. @§. 



Analogous considerations are valid in the on-shell evaluation of Ar (see Section 5 and 
Ref. The threshold corrections that are not suppressed by reduced-mass effects occur 

again in 11^(0) and Tie A z l(m 2 z )/m 2 , the latter being by far the dominant contribution for 
large m t . Here one sets s 2 — > s 2 and c 2 — > c 2 everywhere. We also recall that in Ar the 
contribution of Eq. ( |32a| ) is enhanced by a factor c 2 /s 2 . 



We end this section with the observation that the magnitude of the IS contribution to 
XmU v (s,m t ,m t )/s can be roughly understood on the basis of a simplified "Bohr-atom" 
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model. We recall that, neglecting hard-gluon corrections, this contribution is approximately 
given by 

ImUV ^ m) ^3 lRl ^ )l \ (s-M^), (33) 

where -Ri,o(0) is the radial wave function at the origin and Mq « 2m t is the mass of the 
IS resonance, 6. In Ref. [p^j , |i?i j0 (0)| 2 and M# have been studied in detail on the basis of 
the Richardson potential. Suppose now that we attempt to estimate this effect using the 
"Bohr-atom" picture. In the case of the one-electron atom, (i? lj0 (0)) 2 = 4(a/z) 3 , where \x is 
the reduced mass. For toponium we set /i = m t /2, Me = 2m t , and replace a — > 4d s (fci)/3, 
where k\ = 2a s (ki)m t /3 is the inverse Bohr radius |2jJ. This leads to 3|_Ri i0 (0)| 2 / (M e m 2 ) = 
(16/9) (a s (&i)) 3 m 2 /m|, which is to be compared with 3m t ^{m t )/ (8m 2 ), obtained from the 
Richardson potential [^4|,|3T|. Here £(mt) is a monotonically increasing function of mt, which 
is evaluated numerically and varies from 1.95 to 3.08 GeV for m z < m t < 250 GeV. For 
m t = m z , 250 GeV, the Richardson potential gives 0.0080 and 0.0347, respectively, while in 
the "Bohr-atom" model the corresponding values are 0.0096 and 0.0428, which are about 
20% larger. This may be partly due to the fact that the Richardson potential is softer 
than Coulombic near the origin. Interestingly, the ratio of the values at m t = 250 GeV and 
m t = m z , which gives an indication of the m t dependence, is almost the same in both cases. 
In the "Bohr-atom" model, for the nS states, |i? n)0 (0)| 2 scales as n~ 3 . However, as the radii 
of the higher orbits increase as n 2 , the momentum k n at which a s is to be evaluated becomes 
smaller. For sufficiently high n, the relevant a s falls in the non-perturbative regime and the 
"Coulombic" picture becomes increasingly doubtful. The evaluation based on a realistic, 
"confining" potential is clearly preferable. However, it is interesting that a simple Coulomb 
potential gives a similar answer, at least in the case of the IS state. 



7 Numerical Calculations 

In the previous sections we have discussed the theoretical background necessary for the 
incorporation of the leading QCD effects in the basic corrections Af w , Af, and Ar. At the 
same time, as explained in Section 1, we have introduced a number of refinements in the 
analysis of the electroweak corrections. 

In this section, we apply the previous results to numerically evaluate Af w , Af, and 
Ar, and, most importantly, to derive precise values for m w and s 2 = sm 2 9 w (m z ), as 
functions of m t and m H . Working first in the ms scheme, the basic strategy is the same 
as in Ref. [Ijj. We employ Eqs. ( |10cj ) and ( P] , |18a| ) as the basic expressions for Af w and 
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Af, respectively, and iteratively evaluate these corrections in conjunction with Eq. ( |19a| ), 
leading to accurate values for s 2 . Then Ar can be evaluated from Eqs. (|19cj) or ( |19e| ), s 2 from 
Eq. (pb|) , and 77i w from either s 2 or Eq. ( |191| ) . We use as input values ol — (137.0359895) 
Gf, = 1.16639 x 10- 5 GeV~ 2 g|, m z = 91.187 GeV @, leading to A = (vm/v^G^) 1 / 2 = 



37.2802 GeV, and d s (m z ) = 0.118 [39]. As shown in Appendix C, the effects arising from 
finite fermion masses are very minor. Nonetheless, we include them as follows: the u, d, and 
s quarks are treated as massless, while we employ m c = 1.5 GeV and mb = 4.5 GeV. The 



leptons are given their physical masses p8[ , including the recent value m T = 1.777 GeV [plfl . 
In order to incorporate the new results on the leading irreducible corrections of 0(a 2 ) [^], 
it is very convenient to use a precise analytical representation of the function R(m H /m t ) in 
Eqs. ( |17a| , |18a| . [23b| ). In the range r = m H /m t > 4, the authors of Ref. J7| give the accurate 



asymptotic expansion 
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R(r) = — + 7r 2 -271nr + 61n 2 r 
4 

+At(2 - 12tt 2 - 24 In r - 108 In 2 r) 

3 r 2 

' (1613 - 240vr 2 + 3000 lnr - 2880 In 2 r). (34a) 
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In the complementary domain < r < 4, only numerical values are available f7j, which we 
have fitted with the expression 

R( r ) = -0.7392088 + r(-11.5315 + 0.3824971nr) 

+ r 2 (5.31338 - 3.055 lnr + 0.523039 In 2 r). (34b) 



The maximum deviation of Eq. (|34b|) from the original data is < 0.025 and occurs near 
the matching point, r = 4. 

In Tables I-V we display the calculated values of m w and s 2 , as functions of m t for 
m H = (60, 100, 250, 600, 1000) GeV. Currently, a fit to all data, using the electroweak radia- 



tive corrections of the SM and a s {m z ) = 0.120±0.006, gives m t = 145+^:^ GeV 0. The 
central value corresponds to m H = 300 GeV and the last error reflects the theoretical uncer- 
tainty associated with the range 60 GeV < m H < 1 TeV. This implies m t < 173 GeV at 
the 95% confidence level, where the error is again due to the m H uncertainty. Although this 
strongly suggests that m t ^ 200 GeV, in the tables we present values up to m t = 250 GeV. 
One of the reasons is that it is theoretically interesting to study and compare the high-m 4 
dependence of the various radiative corrections. Moreover, there is always the possibility 
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of a statistical surprise or that unknown new physics may alter the predictions of the SM, 
so that scenarios in which m t is found at higher values are not completely excluded. For 
comparison purposes, we list in the tables the results of four different calculations: i) only 
electroweak corrections, with all QCD corrections turned off (columns labelled EW for 
"electroweak" ) ; it) electroweak plus perturbative (D(aa s ) corrections (columns labelled P 
for "perturbative"); Hi) the above, plus threshold effects calculated in the resonance ap- 
proach according to Section 6 (columns labelled P+R for "perturbative plus resonance"); 
iv) the same as in (iii) but with threshold effects evaluated in the Green-function approach 
(columns labelled P+G for "perturbative plus Green function"). In order to keep the tables 
compact, we have displayed only the quantities of greatest physical interest, namely m w 
and s 2 , rather than the radiative corrections Af w , Af, and Ar, or the derived parameter s 2 
of the on-shell scheme. The interested reader can readily glean these important quantities 
from the tables. Thus, inserting the value of s 2 in Eq. and those of s 2 and m w in 
Eq. (^), one finds Af and Af w , respectively. Using m w , one calculates s 2 = 1 — m^/m 2 
and, in conjunction with Eq. ([[]), Ar. We have kept enough decimal figures in s 2 and 
so that Af w , Af, Ar, and s 2 can be accurately evaluated. The tables allow us to 



w 1 



separate the threshold effects from the more established perturbative 0(aa s ) contribu- 
tions. Although the resonance (R) and Green-function (G) approaches are quite different 
conceptually and technically, the tables reveal the welcome and rather surprising result 
that their effect on m w and s 2 is very similar over the entire range m z < m t < 250 GeV, 
60 GeV < m H < 1 TeV. Nonetheless, for reasons explained in Ref. [31| and Section 6, in 
our specific applications we use the resonance method for m t < 130 GeV and the Green- 
function approach for m t > 130 GeV. 

It is also a curious and rather surprising fact that most of the new effects we have 
considered relative to Ref. |1[, namely a) the incorporation of the recent results on the 
leading irreducible corrections of 0(a 2 ), b) the perturbative 0(aa s ) corrections, and c) the 
threshold contributions in the formulation of Section 6, increase the values of Af and Ar 
for given m t and m H . Thus, they have a sign opposite to that of the leading m r dependent 
part of the one-loop corrections and, therefore, they induce an increase in the mt upper 
bound. For large mt, they are only partially compensated by the shift of —5 x 10 -4 in 
Af and Ar arising from the new calculation of e 2 7Ze (n( 5 7 )(o)-n( 5 7 )(m|)) |pj. it should 
also be remembered that most of the corresponding shifts in Af and Ar increase with 
m t : (a) behaves as m\ (cf. Eqs. (|17a| , |18a| , |23b| )), the dominant (b) contributions as m 2 (cf. 
Eqs. ( |16a| , |19c|) ), and (c) more rapidly than linearly in m t ||31|| . 

The QCD effects on m w and s 2 are visible in the tables. For example, for the interme- 
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diate value m H = 250 GeV and m t = m z , 130, 150, 180, 200, 250 GeV, the perturbative 
0(aa s ) corrections lead to the shifts 



Am w = -(32,42,50,65,77, 113) MeV 
As 2 = +(1.4,2.3,2.8,3.7,4.4,6.4) x 10~ 4 

The threshold contributions lead to further shifts 

Am w = -(8,14,14,20,25,42) MeV 
As 2 = +(0.7, 0.9, 0.9, 1.2, 1.5, 2.5) x 10~ 4 



Q{aa s ). (35a) 



(threshold). (35b) 



The departure from the monotonic behaviour between mt = 130 GeV and m t = 150 GeV is 
due to the transition from the resonance to the Green-function approach. The effect of the 
leading irreducible corrections of 0((a mf/m^) 2 ) cannot be seen in the tables because they 
have been included in every column. However, they can be estimated from the relations 

Am w _ c 2 R (H?2L \ ^± 



m w c 2 — s 2 — 2c 2 x t V m t J 6 ' ^ ^ 

A* 2 - -4^, R (—) W^U' (36b) 
c — s V m t J 3(1 — x t ) 



where Xt is defined in Eq. ( |17b|) . The term 2c 2 Xt in the first denominator of Eq. ( |36a| ) 
and the factor (1 — x t )~ 2 in Eq. ( |36b| ) take into account the fact that some of the leading 
one- loop contributions to Ar and Af depend on s 2 and s 2 , respectively, and are therefore 
affected by the shifts in these parameters. For the same values of m H and m t as employed 
before, Eqs. ( |36a| ,b) give 



Am w = -(1,5,8, 16,24,52) MeV 
As 2 = +(0.1,0.3,0.5,0.9,1.4,3.0) x 10~ 4 



0(a 2 ). (36c) 



Equations ( p5a| ,b, ^6^ ) can be compared with the experimental uncertainties (8 



m 



w ) cxp 



100 MeV and (5s 2 ) exp w 4 x 10" 4 expected at the end of 1993 HJ. We see that, in general 



terms, they are of the same order of magnitude. Of course, in the long run, even better 
experimental accuracies are envisaged, reaching perhaps (Sm w ) exp w 50 MeV. The above 
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shifts can also be compared with the theoretical uncertainties 

(fav), t = g2 _ /_ 2e%< 1 « ±16 MeV, (37a) 



c 2 s 2 £Af 



(^ 2 ) th = (1^2) j^a; w ±3 x 10 " 4 ' (37b) 



arising; from the error (<5Ar) th « (5Ar) th = ±9 x 10~ 4 in the calculation of e 2 7& (n<g(0) 
— n^(m|)J [pTO | . We recall that the theoretical error arising from the neglect of higher- 



order electroweak corrections is expected to be of 0((a/irs 2 )xt) in Af w and Af and further 
enhanced by a factor c 2 /s 2 in Ar (see Section 1). Moreover, the error in the overall 



QCD corrections (perturbative and threshold effects) has been estimated to be ±20% 
Incidentally, Eq. (|37b| ) shows that, if the experimental accuracy in s 2 is improved in the 
future well beyond 4 x 10~ 4 , a meaningful theoretical interpretation will require a decrease 
in the above-mentioned theoretical errors. 

As pointed out before, the higher-order corrections we have considered lead to an in- 
crease in the mt values obtained from experiments. As an example, we consider the case 
m H = 250 GeV and call m t the parameter derived from (m w ) cxp when the perturbative 
0(aa s ) corrections are included. Table III shows that the possible values m t = m z , 130, 
150, 180, 200 GeV are larger than those obtained in the purely EW calculation by 

Amt = +(5.8, 7.5, 8.3, 9.7, 10.8) GeV (0{aa s )). (38a) 

The values of mt derived from (s 2 ) ex p are shifted by slightly higher amounts, the differences 
with Eq.( [38a|) being <0.8 GeV. The bulk of Am t arising from the perturbative 0(aa s ) 
corrections can be understood with the simple formula ( |16c| ), which describes the dominant 
contribution. Similarly, we see from Table III that, when threshold contributions are in- 
cluded, the values of m t derived from (m w ) exp are larger than those obtained in the EW+P 
calculation by additional shifts of 

Am t = +(1.7, 2.7, 2.7, 3.3, 4.1) GeV (threshold). (38b) 

The corresponding variations arising from the irreducible 0(a 2 ) corrections are 

Am t = +(0.6, 1.2, 1.6, 2.5, 3.3) GeV (0(a 2 )). (38c) 

Thus, for m H = 250 GeV and m t = 200 GeV, the combination of (D(aa s ), threshold, and 
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leading irreducible 0(a 2 ) corrections increases the value of m t derived from (m w ) cxp by 
« 16.8 GeV. 

Another topic of considerable interest is the comparison of calculations carried out in the 
Is and on-shell methods of renormalization. This is illustrated for m H = 60, 250, 1000 GeV 
in Tables VI- VIII, where we show the evaluation of Ar, a physical observable, obtained 
on the basis of the Is approach of Section 3 and the on-shell formulation of Section 5. 
For the purposes of this study, we have included the electroweak and perturbative 0(aa s ) 
corrections, leaving aside the threshold effects. We also display the derived values of m w and 
s 2 . Inspection of the tables shows that the two calculations of Ar are in excellent agreement 
over the entire range 60 GeV < m H < 1 TeV, m z < m t < 250 GeV, with a maximum 
discrepancy of ~ 2.5 x 10~ 4 occurring at m t = 250 GeV and m H = 1 TeV. Although such 
accurate agreement may be somewhat fortuitous, it is roughly of the expected order of 
magnitude, i.e. 0((a/rrs 2 )(c 2 /s 2 )x t ). 
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Appendix A 



In this Appendix we discuss the perturbative corrections of 0(a s ) to the vacuum polar- 
ization functions involving quarks. Defining U v,A (s, m 1 , m 2 ) according to Eqs. ( |26a| ,b), we 
expand 

(A.l) 



n 



V,A< 



s, mi,m 2 ) 



n, 



V.A, 



s, mi, m 2 ) H Hi ( s > m ii m 2 
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where s = q 2 and mi and m 2 are the masses of the virtual quarks in the loop. The functions 
Yl^' A (s, mi, m 2 ) have been extensively discussed in the literature. They can be gleaned, for 
example, from Ref. fl4~I |. In the 0(a s ) terms we consider two limiting cases: mi = m 2 = m 
(as occurs in the 77, ZZ, and Z7 self energies) and mi = m, m 2 = (as applies, to a very 
good approximation, to the (t, b) contribution to the WW self energy). 
Comparison of Refs. [^2| and [^] leads to the following expressions: 
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(A.2) 



(A.3) 
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where r = s/(4m 2 ), x = s/m 2 , I = ln(m 2 /^') (// is the rescaled 't Hooft mass discussed 
in Section 2), Vi(r), Ai{r), and F\{x) are complicated functions studied in Ref. |2~5| , and 
the color factor appropriate to iV c = 3 has been explicitly included. We recall that ((2) = 
7r 2 /6 and ((3) = 1.20206. . . . The above expressions can be used to evaluate the 0(a s ) 
corrections employed in the text. Thus, the quantity 11^(0) in Section 2 can be written as 



n#(o) = n«(o) + n«(o) + -Reu^(m 2 z ) + -Re (n$(o) - 



77 



77 



~ 77 (™1 



(A.5) 



where the superscripts (/), (t), and (5) refer to the contributions of the leptons, the top 
quark, and the first five quark flavors, respectively. It is easy to see that the contribution 
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Ilj$(s) of quark q to n 77 (s) is 



Q^Ir (s, m q ,m q ) 



(A.6) 



where Q g is the charge of the quark in units of the positron charge e, and m q is its mass. 
In particular, 



nW(0) = Q?(n v ) / (0,m tl m t ) 



(A.7) 



where the prime on the r.h.s. denotes differentiation with respect to s. Using our Eq. (|A.2| ) 
and Eq. (18) of Ref. [[25J, we find that the 0{a a ) part of II$(0) is (2d s /9vr 3 ) [ln(/i'/m t ) 
+ 15/8 — (2(n — 4)) -1 ]. The finite part of this result is included in the second term of 
Eq. (|7aD , while the pole contributes to the last. Similarly, 



7&n$(m|; 



ST* 



(A. 



Using Eq. (TAg ), we see that the £>(a s ) part of Eq. flAlp is (d s /4vr 3 ) E^Qj [ 2 ln(/x'/ 
+f2(r q ) - ( 



(n — 4) _1 ], where r q and /a(r) are defined in Section 2. The finite and pole parts 
have been included in the third and last terms of Eq. (|7a|), respectively. As mentioned in 
Section 2, for e 2 7Ze ^11^(0) — II^(m|)j we employ a recent evaluation [0]. The contri- 
butions of 0(a) within the square brackets in the second and third terms of Eq. ( [7a] ) are 
simply obtained from the 0(a s ) ones by dividing by the quadratic Casimir coefficient 4/3 
for the fundamental representation of SU(3)c, multiplying by an additional factor of Q 2 q1 
and replacing a s — > a. The contributions of 0(a) in 11^,(0) (first term in Eq. fl7a|)) can 
be obtained from those in Iiw(0) by changing m t — > m ; and dividing by 3Qf , where the 3 
stands for the color factor. 

In the approximation of neglecting the squares of mixing angles, the quark contribution 

tO Ay/w IS 
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In the 0(a s ) part we approximate mi = ui2 = in the (u, d) and (c, s) contributions and 
mi = mt, = in the (t, b) contribution. In this limit, H A (s, mi, 7712) = Il y (s,mi,m2) 
and we see that the 0(a s ) part of Eq. ( |A.8| ) is 
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where the first and the second terms correspond to the contributions of the light and (t, b) 
isodoublets, respectively. The choice of renormalization scale in a s was explained at the 
end of Section 2. Employing our Eq. QA.4Q and Eq. (20) of Ref. [p5l, we find 
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Subtracting the pole term, setting // = m z , and multiplying by e 2 /s 2 , we obtain Eq. ( |lld|) . 
Similarly, the quark contribution to A zz (s) is 
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where C 3g = +1 (—1) for up (down) members of the doublet and the sum is over quarks. 
In the 0(a s ) part we neglect all masses other than m t and find 
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where the first and second terms are the contributions of the light and (t, b) isodoublets, 



35 



respectively. Combining Eqs. ( |A.1(J| , 



, A.14| ) and recalling Eqs. ( |lld|jA.2| -4), we obtain 
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Subtracting again the pole terms, setting // = m z , and multiplying by e 2 /s 2 , we obtain 
Eq. (1153) . 

Finally, in order to implement the decoupling of the top quark, we need the quantity 
(a/n)d, where d is defined in Eq. (|9l5D. From Ref. |41[] we find 
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3 V2s 2 3, 

The 0(a s ) part of Eq. ( |A.17| ) is evaluated as in the case of 11^(0) (cf. discussion after 
Eq. ( |A.7| )). A relevant combination that occurs in Af w and Af is 
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where d is defined in Eq. (I0d). 



Appendix B 

In this appendix we outline the derivation of the expressions for Af w (cf. Eq. (|10a|) ) and 
Af (cf. Eq. (|I^)), corresponding to the strategy in which the decoupling of the top quark 
is implemented. 

As s 2 (l — Af w ) and s 2 c 2 (l — Af) are physical observables (cf. Eqs. (@|f)), we define 
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Af w and Af in such a way that 



[S 2 (l-Af w )] old = [S 2 (l-Af w -)] ncw , 
[^ 2 (l-Af)] old = [S 2 c 2 (l-Af)] ncw . 



(B.l) 
(B.2) 



Here the subscript "old" labels the quantities obtained when the Is counterterms cancel 
only the divergent parts involving 5 (cf. Eq. Q5ap), as in Ref. M, while new denotes the 
corrections employed in the present paper, where the Ws counterterms contain small fi- 
nite parts necessary to implement the top-quark decoupling in relevant amplitudes. Our 
strategy is to retain, in Eqs. ( B.l ) and ( |B.2| ), terms of 0(a 2 ) when they involve large log- 
arithmic or mf/m 2 enhancements, but to neglect them otherwise. We also neglect certain 



corrections of 0(a 3 ). Recalling = e 2 /Z e and Eqs. (|7d| , pa| ,b, |A.18|) , we have 
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where d and d are defined in Eqs. ( Pc|) and ( |10d| ), respectively. 

The derivation of Eq. ( |10a| ) follows almost immediately. Inserting Eq. ( p.3| ) into 
Eq. (pll) , we find 
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As the only large correction in Af w is (a/7r)A 7 and d(l — aA 7 /7r) = a (cf. Eq. flSc])), we 
obtain 



(Af v 
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(B.6) 



The expression for (Ar w ) Q id in terms of 11^(0), WW self energies, and the vertex- and box- 
diagram corrections to \x decay, is given in Eqs. (7b, 8b) of Ref. The only contribution 
to (Arv)oid with large logarithmic enhancement is e 2 11^(0) \jfs and, to 0(a), this does 
not involve s 2 . Therefore, we can replace s 2 ld — > s 2 ew everywhere in (Af w ) old , neglecting 
very small terms of 0(a 2 ). Inserting the expression for (Af w ) id m Eq. ( p.6| ), we obtain 
Eq. 0. 



The derivation of Eq. (14) is more subtle. Inserting Eq. ( |B.3|) and 

d] , 
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into Eq. ( |B.2| ), we have 
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We first consider the second term in Eq. ( |B.8[ ). Using Eq. ( |18b| ), this contribution can 
be written as —{a/^)d{l — [s 2 /c 2 ] new ) [1 — Af w ] old c 2 /c 2 ld , which, neglecting small 0(a 2 ) 
terms, becomes — (a/ir)d (1 — [S 2 /c 2 ] new ) c 2 /c new . Next we analyse the first term on the 
r.h.s. of Eq. (|B.8|). According to Eq. (15b) of Ref. D, it is given by 
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Examination of Eq. (A. 8) of Ref. shows that 
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where the ellipses represent non-leading terms (some of which involve ln(mt/m z ) with very 
small coefficients proportional to s 2 , s 2 , or s 4 ). The significant point is that the leading 
contributions in Eq. ( [B . 1 0|) are independent of s 2 . Therefore, in Eq. (p.9| ) we can replace the 
last factor in the second term by an analogous expression with s 2 ld — > s 2 ew , 



s 2 ,„, the difference 



being again small terms of 0(a 2 ). Using Eqs. ( |B.4|JB.6| , [B.7| ), Eq. ( |B.9| ) becomes 
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0(a 2 ,a 3 ), we find that the term proportional to (a/n)d(l — s 2 /c 2 ) in Eq. ( |B.11| ) can be 
written as (a/ir)d (1 — s 2 /c 2 ) (c 2 /c 2 — 1). Inserting the above results into Eq. ( |B.8|) , we find 
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The coefficient of (a/7i)d simplifies to s 2 /c 2 and Eq. ( |B.12|) equals Eq. (0). 
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Appendix C 



If terms proportional to squares of mixing angles are neglected, the fermionic contributions 
to the WW self energy become a sum over independent isodoublet contributions and we 
obtain 
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where B^ is defined in Eqs. ( |lla| ,b), N c = 3 and N c — 1 for quark and lepton isodoublets, 
respectively, and cj-t- = vn? ± jrr? wl with m + and m_ being the masses of the "up" and "down" 
fermions in the isodoublet. Calling 
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the function Q(u + , uJ) is given by 
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Equations ( G.4 , |U.5|) are a simpler version of Eqs. (A. 10, A. 11) of Ref. [fi9 ]. As expected 
from the integral representation of the self energies JIT], fi(ui + , is a symmetric function 
of u; + and u;_. As c<j_ — > 0, 



/ 1 (^ + ,0) = ln(cV) + (^ + -l) 2 
while in the limit uj + ,uj_ — > 0, 
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If all the fermion masses other than m t are neglected, Eqs. (|C.1 , C.6 , |C.7 ) lead to Eq. (|llcj) . 
In order to estimate the terms of 0(myTn^), we consider the case uo + > 1 and uj- <C 1. 
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Including corrections of 0(u)-), we have 
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Inserting the second term of Eq. ( |C.8|) in the (t, fe)-isodoublet contribution to Eq. ( |C.1|) 
and identifying oj- = uj b = mj/mj, and uj + = uj t = m^/m^, we obtain a very small shift, 
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Using m b = 4.5 GeV and m z = 91.187 GeV, we find that the 0(ml/m\,) corrections to 
Eq. (|CTTD are « 2.5 x 10~ 5 , 5.4 x 10" 6 , 2.8 x 1(T 6 , 1.7 x 10~ 6 for m t = m z , 150, 200, 
250 GeV, respectively. We note that the cofactor of u b in Eq. ( p.9|) tends to as m t — > oo. 
To estimate the contributions of C^m^/m^), we set m s = and keep terms of 0(uj + ) for 
u + « 1 in Eq. Q: 
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Inserting the second term of Eq. (|C.10|) in the (c, s)-isodoublet contribution to Eq. ( |G1|) , 
and identifying u + = uj c = rv? c jrr? w ^ we obtain 
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For m c = 1.5 GeV, this amounts to a shift of —1.1 x 10~ 5 to Eq. ( |C.1| ). For the corrections 
of 0{m 2 ./m^ v ), we replace uj c — > co> T = m^/m^, and divide by the color factor 3 in Eq. ( |C11| ). 
This leads to a further correction of —5 x 10~ 6 to Eq. ( |G1|) . 

Turning our attention to Cq , defined in Eqs. (|15a| ,b) , and neglecting again the squares 
of mixing angles, we have 



(/) 



a 



12ns 2 c 2 



c 2 J2 N c f 2 (u + ,U-)-W N c9(»f) 

doublets / 



(C.12) 



40 



1 - 2^+ + ^-) 
x2fl(w + , u_) 



1 - 



(Jj-L. — U>- 



In + uj + + u- + 



(w + — 



g(ti f ) = [(1 - 4s 2 C 3f Q f ) 2 + l] [in/x, + - + 2(1 + 2 N ) (A(D f ) - 1) 
-6/// (ln// / + 2A( J D / )-2), 



(C.13) 



(C.14) 



where the / sum is over quark and lepton flavors, N c is again the color factor, Qf is the 
electric charge in units of the positron charge e, C 3 f — +1 (— 1) for up (down) members of 
the doublet, /i/ = m 2 /m 2 , Df = 4/i/ — 1, and 
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As is clear from their structure, the contributions involving f 2 and g in Eq. flC.12| ) arise 
from the first and second terms in Eq. ( |1 5a| ) , respectively. 
One readily finds the limiting values 
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If the fermion masses other that m t are neglected, Eqs. ( |C . 1 2| , PTT4| , |C . 17] — 19) lead to Eq. (|15cj ). 
To discuss the corrections of O^ml/m^), we consider / 2 (a>_|_, a>_) for w + > 1, w_ < 1 and 
<?(///) for ii f -C 1. One readily finds 
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= s(0)- 4&n f &Q f {2&Qi -C 3/ ) + C(// 2 / ). (C.21) 
Identifying u;_ = cj&, u + = u t , [if = fit, = ml/m 2 and inserting Eqs. (|CJ . 2 0| , |CT72T| ) into 
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Eq. ( |(J.12p , we find for the leading correction of 
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Using the same input values as in Eq. (|C.9|) , we find that the 0{m 2 /m 2 v ) corrections to 
Eq. flCTT^ ) are « 2.5 x 10~ 5 , -6.7 x 10~ 6 , -1.6 x KT 5 , -2.3 x 10" 5 for m t = m z , 150, 200, 
250 GeV, respectively. To estimate the corrections of 0{m 2 /m 2 v ), we set u_ = and keep 
terms of 0(u+) for u + <C 1 in Eq. ( |C. 17|) . In conjunction with Eq. ( |C.21|) , this gives 
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where fi c = m 2 c /m? z . For m c = 1.5 GeV, this is ~ —8 x 10 -7 . The corrections of 0(m 2 ./m 2 v ) 
are even smaller on account of the absence of the color factor. 

There are also very small corrections of 0(m 2 /m 2 z ) and 0(m 2 /m 2 ) associated with the 
fi{r q ) contribution to e 2 IT^(0) (cf. Eq. ©). Their magnitude is ~ +5 x 10 . 

Putting all these small 0(jn 2 /m 2 v s j corrections together, we see that for m t = m z , 150, 
200, 250 GeV they amount to (1.4, -0.6, -0.8, -0.9) x 10~ 5 in the case of Af w and (-1.0, 
+0.2, +0.9, +1.5) x 10~ 5 in the case of Af. On the other hand, Ar = Af w — (c 2 /s 2 )(Af w — 
Af) and the small corrections are enhanced in the second term, leading to shifts of (—6.5, 
+2.1, +5.1, +7.8) x 10^ 5 for the same values of m t . 
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22853 


81 
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80 


916 
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80 


878 


250 





22721 





22785 


0.22806 





22809 


81 


114 


81 


001 


80.964 


80 


959 



Table I: Calculated values of sin 6 w (m z ) and Til w , clS cl function of m t , for m z = 
91.187 GeV, a s (m z ) = 0.118, and m H = 60 GeV. The EW column includes only elec- 
troweak radiative corrections. The P column incorporates perturbative 0(aa s ) contribu- 
tions (see Section 3). The P+R and P+G columns contain also ti threshold effects (see 
Section 6) in the resonance and Green-function approaches, respectively. The on-shell pa- 
rameter sm 2 9 W and the radiative corrections Af w , Af, and Ar can be gleaned from this 
table, as explained in Section 7. 
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Table II: As Table I, for m H = 100 GeV. 
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Table III: As Table I, for m H = 250 GeV. 
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80 
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80 
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80 
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80 
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80 
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80 


803 


80.765 


80 
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Table IV: As Table I, for m H = 600 GeV. 
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79 
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23329 
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150 
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23326 


0.23337 
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160 
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80 
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80 
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23193 





23230 
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80 
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259 
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80 


239 


190 





23155 





23196 


0.23211 





23209 


80 


395 


80 


324 


80.299 


80 


302 


200 





23116 





23160 


0.23176 





23175 


80 


468 


80 


390 


80.363 


80 


366 


210 





23075 





23123 


0.23141 





23140 


80 


543 


80 


459 


80.430 


80 


431 


220 





23033 





23085 


0.23104 





23104 


80 


621 


80 


530 


80.499 


80 


499 


230 





22990 





23046 


0.23066 





23066 


80 


700 


80 


602 


80.569 


80 


568 


240 





22946 





23006 


0.23027 





23028 


80 


782 


80 


677 


80.641 


80 


638 


250 





22900 





22965 


0.22987 





22990 


80 


866 


80 


753 


80.715 


80 


710 



Table V: As Table I, for m H = 1000 GeV. 
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[GeV] 




m z 


6.003 


6.006 


79.947 


79 


946 


100 


5.735 


5.738 


79.996 


79 


995 


110 


5.437 


5.440 


80.050 


80 


049 


120 


5.135 


5.138 


80.103 


80 


103 


130 


4.823 


4.826 


80.158 


80 


158 


i a n 
14U 


4.498 


4.502 


80.215 


80 


O 1 A 
2,14: 


150 


4.159 


4.162 


80.274 


80 


273 


160 


3.801 


3.805 


80.334 


80 


334 


170 


3.425 


3.429 


80.398 


80 


397 


180 


3.029 


3.033 


80.463 


80 


463 


190 


2.611 


2.615 


80.532 


80 


531 


200 


2.170 


2.173 


80.603 


80 


602 


210 


1.705 


1.708 


80.677 


80 


676 


220 


1.214 


1.217 


80.754 


80 


753 


230 


0.697 


0.699 


80.833 


80 


833 


240 


0.151 


0.153 


80.916 


80 


915 


250 


-0.424 


-0.423 


81.001 


81 


001 



Table VI: Comparison between the values of Ar and m w obtained using the ~ms approach 
of Section 3 ((Ar)/, (m w )j) and the on-shell formulation of Section 5 ((Ar)//, (m w )jj), as 
a function of m t , for m z = 91.187 GeV and m H = 60 GeV. Nonperturbative tt threshold 
effects are not included here. 
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(m 
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(m 


w)ll 


[GeV] 
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[GeV] 




m z 


6.461 


6.462 


79 


862 


79 


862 


100 


6.199 


6.200 


79 


911 


79 


911 


110 


5.907 


5.908 


79 


965 


79 


964 


120 


5.612 


5.613 


80 


018 


80 


018 


130 


5.309 


5.309 


80 


072 


80 


072 


140 


4.994 


4.994 


80 


128 


80 


128 


150 


4.665 


4.665 


80 


186 


80 


186 


160 


4.320 


4.320 


80 


246 


80 


246 


170 


3.959 


3.957 


80 


308 


80 


308 


180 


3.579 


3.577 


80 


372 


80 


372 


190 


3.179 


3.176 


80 


439 


80 


439 


200 


2.760 


2.756 


80 


508 


80 


508 


210 


2.319 


2.313 


80 


579 


80 


580 


220 


1.855 


1.848 


80 


653 


80 


654 


230 


1.369 


1.360 


80 


730 


80 


731 


240 


0.858 


0.846 


80 


809 


80 


810 


250 


0.322 


0.307 


80 


890 


80 


892 



Table VII: As Table VI, for m H = 250 GeV. 
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7.036 


7.037 


79 


754 


79 


754 


100 


6.779 


6.780 


79 


803 


79 


803 


110 


6.494 


6.494 


79 


856 


79 


856 


120 


6.206 


6.205 


79 


910 


79 


910 


130 


5.910 


5.910 


79 


964 


79 


964 


140 


5.604 


5.603 


80 


019 


80 


020 


150 


5.286 


5.284 


80 


077 


80 


077 


160 


4.954 


4.951 


80 


135 


80 


136 


170 


4.607 


4.603 


80 


196 


80 


197 


180 


4.244 


4.239 


80 


259 


80 


260 


190 


3.865 


3.859 


80 


324 


80 


325 


200 


3.469 


3.461 


80 


390 


80 


392 


210 


3.056 


3.045 


80 


459 


80 


461 


220 


2.625 


2.612 


80 


530 


80 


532 


230 


2.175 


2.159 


80 


602 


80 


605 


240 


1.707 


1.687 


80 


677 


80 


680 


250 


1.220 


1.195 


80 


753 


80 


757 



Table VIII: As Table VI, for m H = 1000 GeV. 
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